Constraining the neutron-proton effective mass splitting in neutron-rich matter 
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Within Bombaci's phenomenological single-nucleon potential model we study the neutron-proton 
effective mass splitting m* — m* in neutron-rich matter. It is shown that an effective mass splitting 
of m* < nip leads to a symmetry potential that is inconsistent with the energy dependence of the 
Lane potential constrained by the nucleon-nucleus scattering experimental data. 
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One of the most fundamental properties characterizing 
the propagation of a nucleon in nuclear medium is its ef- 
fective mass 0,0, El- The latter describes to leading order 
the effects related to the non-locality of the underlying 
nuclear effective interactions and the Pauli exchange ef- 
fects in many-fermion systems. In neutron-rich matter, 
an interesting new question arises as to whether the ef- 
fective mass m* for neutrons is higher or lower than that 
for protons to*. Microscopic many-body theories, e.g., 
the Landau-Fermi liquid theory|4| and the Brueckner- 
Hartree-Fock (BHF) approach [3, predict that to* > to* 
in neutron-rich matter. On the other hand, some effective 
interactions within phenomenological approaches predict 
the opposite 6]. Unfortunately, up to now almost noth- 
ing is known experimentally about the neutron-proton 
effective mass splitting to* — to* in neutron-rich medium. 
Knowledge about nucleon effective mass in neutron-rich 
matter is critically important for understanding several 
properties of neutron stars0, H, @- It is also important 
for the reaction dynamics of nuclear collisions induced 
by radioactive nuclei, such as, the degree and speed of 
isospin diffusion, the neutron-proton differential collec- 
tive flow as well as the isospin equilibrium!! G3, EJ E3 
Moreover, it influences the magnitude of shell effects and 
basic properties of nuclei fare from stability [l| . How- 
ever, even the sign of the neutron-proton effective mass 
splitting remains rather controversial theoretically. It is 
thus imperative to constrain the neutron-proton effective 
mass splitting in neutron-rich matter using experimental 
data. In a recent study by Rizzo et al.0, nuclear reac- 
tions with radioactive beams were proposed as a means to 
disentangle the sign of the neutron-proton effective mass 
splitting in neutron-rich matter. In this note we take a 
completely different approach for the same purpose. We 
show that an effective mass splitting of to* < to* leads to 
a symmetry potential that is inconsistent with the energy 
dependence of the Lane potential constrained by existing 
nucleon-nucleus scattering data. The effective mass split- 
ting of to* < to* is thus ruled out phenomenologically 



using experimental data indirectly. 

The effective mass to* of a nucleon r (n or p) is deter- 
mined by the momentum-dependent single nucleon po- 
tential U T via 3] 
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where k^ is the nucleon Fermi wave number in asymmet- 
ric nuclear matter of isospin asymmetry 6 = (p n — Pp)/p 
with p, p n and p p being the nucleon, neutron and proton 
density, respectively. The effective mass defined in Eq. Q] 
is the so-called k-mass which is identical to the total ef- 
fective mass when the single nucleon potential is energy- 
independent 0, 0]. For the single nucleon potential U T , 
following Rizzo et al.|6j we use the phenomenological for- 
malism of Bombaci[l3j 

U r (k,u,5) = Au + Bu a 
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where u = p/po is the reduced density and ± is for neu- 
trons/protons. In the above T T — j^jw J d 3 kf T (k)g(k), 
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g(k) is a momentum regulator g(k) = 1/ 

and f T (k) is the phase space distribution function. The 
parameter A = 1.5Kp where Kp is the nucleon Fermi 
wave number in symmetric nuclear matter at normal 
density p . With A=-144 MeV, B=203.3 MeV, C=-75 
MeV and a = 7/6 the Bombaci formalism reproduces 
all ground state properties including an incompressibil- 
ity of ^0=210 MeV of symmetric nuclear matte[(| fl3j . 
It should be noticed that the Bombaci expression is an 
extension of the well known Gale-Bertsch-Das Gupta 
formalism^! from symmetric to asymmetric nuclear 
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matter. The various terms are motivated b y th e HF anal- 
ysis using the Gogny effective interactional 0] . This 
potential depends explicitly on the momentum but not 
the total energy of nucleons. Thus the k-mass obtained 
is the same as the total effective mass. One should also 
mention that only the last term is momentum dependent. 
Thus the A parameter sets the scale of the momentum 
dependence of the nucleon potential U T and also the scale 
of the effective mass. The value of A parameter was de- 
termined by the ground state properties of symmetric 
nuclear matter. It can also appropriately reproduce the 
momentum dependence of the empirical isoscalar nuclear 
optical potential as we shall show in the followi ng. The 
expressions in Eq. pleads to an effective mass 0,0] 
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It is noticed that the (1 ± 5) 2 / 3 it 2 / 3 term comes from the 
nuclearn Fermi wave number k^ squared, and ± is for 
t = n/p. The three parameters xq,x 3 and z\ can be 
adjusted to mimic different predictions on the density- 
dependent symmetry energy and the neutron-proton ef- 
fective mass splitting. Two sets of parameters can be 
chosen to give two opposite nucleon effective mass split- 
tings, but almost the same symmetry energy E sym (p) [|| . 
The parameter set z\ = —36.75 MeV, xo — —1.477 and 
X3 = —1.01 (case 1) leads to to* > to* while the one 
with zi = 50 MeV, x = 1.589 and x 3 = -0.195 (case 2) 
leads to to* < to* at all non-zero densities and isospin 
asymmetries. Shown in Fig. 1 are the nucleon effective 
masses as functions of density and isospin asymmetry in 
both cases. It is seen that the neutron-proton effective 
mass splittings, opposite in signs, increase in magnitudes 
in both cases with the density and isospin asymmetry. 
Thus large magnitides of effective mass splittings can be 
obtained in dense neutron-rich matter. These two sets of 
parameters were shown to lead to rather different results 
for several observables in heavy-ion reactions induced by 
neutron-rich nuclei 0. 

Our analysis is based on the consideration that both 
the isoscalar and isovector parts of the single nucleon po- 
tential must have asymptotic values at po in agreement 
with the real parts of the corresponding nucleon optical 
potentials constrained by nucleon-nucleus scattering ex- 
periments. We first examine the isoscalar potentials with 
the two sets of parameters. The quantity (U n + U p )/2 
is a good approximation of the isosclar potential since 
the 6 2 in cq. 121 is negligibly small. It is now customary 
and a more stringent test of the isoscalar potentials to 
compare them with the variational many-body (VMB) 

predictions by Wiringa[H[llE[l3. In th<3 VMB tllC " 
ory the single nucleon potential was obtained by using a 

realistic Hamiltonian that fits nucleon-nucleon scattering 

data, few-body nuclear binding energies and nuclear mat- 
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FIG. 1: (Color on line) Nucleon effective masses as a func- 
tion of density (upper window) and isospin asymmetry (lower 
window) with the two parameter sets (see text). 



ter saturation properties. It also reproduces the experi- 
mental nucleon optical potential available mainly at low 
energies|20j. Shown in Fig. 2 is a comparison of isosclar 
potentials using eq. [21 with the VMB predictions using 
the UV1A two-body potential and the UVII three-body 
potential^f|. First of all, the isoscalar potentials are al- 
most independent of the neutron-proton effective mass 
splittings as one expects. In both cases the isoscalar 
potentials using eq. [21 are in good agreement with the 
VMB predictions up to about k — 2.5/to -1 . At higher 
momenta where combinations of different two-body and 
three-body forces lead to somewhat different predictions 
especially at high densities [lfij. the Bombaci formalism 
leads to slightly lower values. Nevertheless, the quality 
of agreement with the VMB predictions found here is 
compatible with those using all other modelsJllGlt. 

We now turn to the isovector part of the nucleon po- 
tential in both cases. The strength of the isovector 
nucleon optical potential, i.e., the symmetry or Lane 
potential|2l|. can be extracted from UL ane = (U n — 
Up)/ 26 at p Q . Systematic analyses of a large number 
of nucleon-nucleus scattering experiments at beam en- 
ergies below about 100 MeV 22] indicates undoubtedly 
that the Lane potential decreases approximately linearly 
with increasing beam energy Ekm, i.e., 



Uhane — a — bE, 



kin 



(4) 



where a ~ 22 - 34 MeV and b ~ 0.1 - 0.2|23j|2J. Shown 
in Fig. 3 are the theoretical symmetry potentials in the 
two cases in comparison with the Lane potential con- 
strained by the experimental data. The vertical bars are 
used to indicate the uncertainties of the coefficients a and 
b in eq. 01 It is seen that with the effective mass splitting 
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FIG. 2: (Color on line) Strengths of the isoscalar potential 
as a function of density at the four wave numbers for case 
1 (lower window) and case 2 (upper window) in comparison 
with the variational many-body calculations 



{ml - m;)/m = ^ + {k p /k n ff7} (5) 



m* > m* (case 1) the strength of the symmetry potential 
decreases with increasing energy. This trend is in agree- 
ment with the experimental indication. Moreover, the 
slope of the calculated symmetry potential with respect 
to energy is about right although the magnitude obtained 
is slightly higher. Since it is not the purpose of this work 
to find the best parameter set to reproduce the empirical 
Lane potential, no attempt is made to readjust any of the 
parameters first proposed in refs.0,0. In case 2, how- 
ever, the most striking feature is that the symmetry po- 
tential increases with the increasing beam energy. This is 
in stark contrast with the experimental indications. The 
characteristically wrong energy dependence of the sym- 
metry potential in this case thus excludes definitely the 
neutron-proton effective mass splitting of to* < m* in 
neutron-rich matter. 

Our finding here is consistent with that of the ear- 
lier work by Sjobergin the framework of the Landau- 
Fermi liquid theoryQ, BHF predictions as well as 
Hartree-Fock calculations using the Gogny effective 
interaction Within the Landau- Fermi liquid theory, 
very interestingly, the nucleon effective mass splitting has 



where /f™, ff p and /™ p are the neutron-neutron, proton- 
proton and neutron-proton quasiparticle interactions 
projected on the I = 1 Legendre polynomial, as for the 
effective mass in a one-component Fermi liquid. Micro- 
scopic NN interactions predict all /i's are negative in 
symmetric matter at nuclear matter density and in asym- 
metric matter at tree-level. It can be seen then that the 
proton effective mass is smaller in neutron-rich matter 
due to the coupling of the protons to the denser neutron 
background, i.e., the term (fc„/fc p ) 2 /" p is dominant in the 
Eq|S] above, leading to m* > to* as shown numerically 
in Fig. 3 of ref.0. 

In summary, using the Bombaci phenomenological for- 
malism for single nucleon potentials in isospin asymmet- 
ric nuclear matter we analyzed the isoscalar and isovector 
nucleon optical potentials in comparison with the em- 
pirical ones constrained by the experimental data. We 
found that a neutron-proton effective mass splitting of 
to* < to* leads to a characteristically wrong energy de- 
pendence of the Lane potential. Thus, the possibility of 
to* < to* in neutron-rich matter is ruled out indirectly 
by the experimental data. 

This research is supported in part by the National 
Science Foundation under grant No. PHY-0088934 and 
PHY-0243571. 



4 



[1] J. Dobaczewski, Acta Phys. Pol. B30, 1647 (1999). 
[2] J.W. Negele and K. Yazaki, Phys. Rev. Lett. 62, 71 
(1981). 

[3] J. P. Jeukenne, A. Lejeune and C. Mahaux, Phys. Rep. 
25, 83 (1976). 

[4] O. Sjoberg, Nucl. Phys. A265, 511 (1976). 

[5] I. Bombaci and U. Lombardo, Phys. Rev. C44, 1892 
(1991); W. Zuo, I. Bombaci and U. Lombardo, Phys. 
Rev. C60, 24605 (1999); W. Zuo, A. Lejeune, U. Lom- 
bardo and J.F. Mathiot, Eur. Phys. J. A14, 469 (2002). 

[6] J. Rizzo, M. Colonna, M. DiToro and V. Greco, Nucl. 
Phys. A732, 202 (2004) and references therein. 

[7] HA. Bethe, Rev. Mod. Phys. 62, 801 (1990). 

[8] J. Cooperstein, Nucl. Phys. A438, 722 (1985). 

[9] M. Farine et al., Nucl. Phys. A696, 396 (2001). 
[10] V.R. Pandharipande and S.C. Pieper, Phys. Rev. C45, 
791 (1992). 

[11] D. Persram and C. Gale, Phys. Rev. C652, 064611 
(2002). 

[12] B.A. Li, C.B Das, Das Gupta and C. Gale, 
nucl-th/0312054 Nucl. Phys. A (2204) in press. 

[13] I. Bombaci, Chap. 2 in Isospin Physics in Heavy-Ion Col- 
lisions at Intermediate Energies, Eds. B. A. Li and W. 
Udo Schroder (Nova Science Publishers, Inc, New York, 
2001). 



[14] C. Gale et al., Phys. Rev. C35, 1666 (1987). 
[15] R.B. Wiringa, Phys. Rev. C38, 2967 (1988). 
[16] C. Gale et al, Phys. Rev. C41, 1545 (1990). 
[17] C.B. Das, S. Das Gupta, C. Gale and B.A. Li, Phys. Rev. 

C67, 034611 (2003). 
[18] P. Danielewicz and G.F. Bertsch, Nucl. Phys. A533, 712 

(1991); Q. Pan and P. Danielewicz, Phys. Rev. Lett. 

70, 2062 (1993); P. Danielewicz, Nucl. Phys. A673, 375 

(2000). 

[19] P. Danielewicz, R. Lacey and W.G. Lynch, Science 298, 
1592 (2002). 

[20] L.P. Csernai, G. Fai, C. Gale and E. Osnes, Phys. Rev. 

C46, 736 (1992). 
[21] A. M. Lane, Nucl. Phys. 35, 676 (1962). 
[22] D.M. Patterson et al., Nucl. Phys. A263, 261 (1976); 

J. Rapaport et al., ibid A313, 1 (1979); ibid, A330, 15 

(1979); R.P. De Vito et al., Phys. Rev. Lett. 47, 628 

(1981); R. De Leo et al., Phys. Lett. B98, 233 (1981); K. 

Kwiatkowski et al., Nucl. Phys. A301, 349 (1978). 
[23] J.R. Rook, Nucl. Phys. A222, 596 (1973); J. Dabrowski 

et al., Can. J. Phys. 52, 1768 (1974); J. P. Jeukenne et 

al., Phys. Rev. CIO, 1391 (1974). 
[24] P.E. Hodgson, The Nucleon Optical Model, 1994 (World 

Scientific) . 

[25] I thank one of the referees for bringing ref.y to my at- 
tention and suggesting the related discussions. 



